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FOREWORD 


This  report  was  prepared  by  Dr.  Frank  W„  Bubb,  Chief  Soientist 
of  the  Office  of  Air  Research,  Hq..  Wright  Air  Development  Center. 
Work  was  completed  under  (Expenditure  Order  Number  I46I-I.  The  re* 
port  is  one  of  a  series  to  be  issued  on  this  projeot.  Others  in 
the  series  will  be  published  as  research  progresses. 
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This  paper  presents  a  new  operational  calculus  applicable  to  the  approximate 
(and  in  the  Halt,  exact)  analysis  and  synthesis  of  linear  physical  systems  (such 
as  servomechanisms,  electrical  olromita,  and  so  on).  Corresponding  to  any  tine 
function  F(t),  a  "^-transform*  F(x)  IB  y(kv)x  is  defined,  P(kv)  being 

•  k*~K 

ordinates  of  F(t)  at  ia&egral  naltipleg  of  a  tiara  interval  v.  The  relation  be¬ 
tween  input  F(t)  to  any  linear  aysftea,  its  output  H(t),  and  its  nemory  function 

(impulsive  adalttance)  H(t),  is  given  by  the  superposition  integral,  H(t)  1= 
rM 

I  P‘<s)H(t-'S)ds.  Corresponding  to  this  convolution  integral,  the  relation  H(x)  — 
-oo 

P(x)H(x)  holds  between  the  respective  P-traneforas.  Conversely,  this  transform 
equation  implies  the  convolution  equation.  This  correspondence  is  identical  lu 
fora  to  the  correspondence  between  convolution  integral  and  its  Laplace  (or 
Fourier)  transform  which  is  ths  basic  theorem  of  the  classical  operational  nethods. 
The  node  of  constructing  a  table  of  P-transforus  is  indicated  by  working  out  the 
transform  of  the  eleeentary  functions  of  aeobanios  and  of  circuit  analysis.  By 
way  of  indicating  how  this  new  operational  calculus  nay  be  developed,  a  set  of 
usoful  transform  theoreas  are  worked  out  for  such  analytical  operations  as  dif¬ 
ferentiation,  integration,  and  so  on.  Finally,  to  exemplify  the  power  of  this 
new  calculus,  a  linear  differential  equation  is  solved,  numerical  results  being 
compared  with  those  of  the  analytical  solution.  An  outstanding  virtue  of  this 
now  operational  calculus  is  that  It  does  not  depart  from  the  time  domain  itself. 

In  a  subsequent  Air  Force  Technical  Report  No.  6086,  this  method  is  used  to  give 
a  simple  explanation  of  Vieuer’a  theory  of  the  smoothing  and  predicting  filter. 
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A  HEW  LINEAR  OPERATIONAL  CALCULUS 


This  paper  presents  ari  cperaticnal  calculus  applicable  to  the  approximate 
(and  in  the  limit,  exact)  analysis  and  synthesis  of  linear  physical  systems, 
such  as  servomechanisms,  electrioal  circuits,  electromechanisms,  and. the  like. 
In  contrast  to  the  methods  based  upon  the  transform  theories  of  Fourier  and 
Laplace,  this  calculus  doeB  not  depart  from  the  time  domain^,  thus  allowing 
a  more  direot  appeal  to  the  intuition. 

The  new  operational  method  takes  origin  in  the  notion  of  the  generating 
function,  used  so  effectively  by  Laplace^- in  the  theory  of  probability. 
Actually,  the  new  method  is  closely  related  to  the  Fourier  and  Laplace  trans¬ 
form  methods,  but  provides  familiar  algebraic  methods  of  executing  analyses 
and  syntheses.  Before  setting  forth  the  new  methods,  we  recapitulate  a  well 
known  bit  of  linear  system  theory  which  brings  out  the  essential  features  of 
any  operational  method. 


THE  SUPERPOSITION  INTEGRAL 

Consider  the  conventional  "black  box"  representation  of  a  physical  sys¬ 
tem  shown  in  Figure  1.  F(t)  represents  the  input  (oause  or  driving  function) 
which  actuates  the  system.  The  output  of  the  system  (effeot  or  response) 
is  represented  by  H(t).  The  system  itself  is  characterised  by  its  MEMORY 
FUNCTION  M(t) „  which  is  the  response  of  the  system  at  time  t  caused  by  a 
unit  impulsive  input  applied  to  the  system  at  time  t  =  0.  M(t)  is  also 
called  the  system  weighting  function  or  its  impulsive  admittance  and  is  the 
Fourier  transform  of  the  system  frequency  response  as  well  as  the* Laplace 
transform  of  the  system  transfer  function.  T-Jnv«-ro<a. 

Input  F(t),  response  H(t)  and  memory  function  M(t)  are  related^  by 
the  superposition  (convolution  or  faltung)  integral 

4 

H(t)  =  (  F(s)M(t-s)ds  (1) 

This  is  easy  to  see  from  the  following  considerations,  illustrated  in  Figure 
2.  The  element  F(s)ds  may  be  regarded  as  a  small  impulse  at  ti»e  s.  Had 
this  been  a  unit  impulse,  the  corresponding  response  at  time  t  (later  by  the 
time  t  -  s)  would  be  M(t-s)  .  The  proportional  response  due  to  the  non-unit 
impulse  is  then  F( s)dsM(t-s) •  The  total  response  H(t)  at  time  t  is  the  sum 
of  all  the  responses  at  time  t  due  to  similar  elementary  impulses  occurring 
preceding  time  t,  namely,  is  the  superposition  integral  (l). 

The  memory  function  M(t)  vanishes  for  negative  t,  since  it  is  axiomatic 
that  an  effeot  or  response  oarrnot  precede  its  cause  (the  unit  impulse  at 
t  =  0) .  This  means  that  M(t-s)  =  0  for  i>  tj  hence  the  upper  limit  t  in 
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(1)  may  changed  to  c*>  .  This  is  also  clear,  since  any  impulse  after  t  cannot 
be  retroactive  to  produoe  any  response  at  the  earlier  time  t. 

It  is  important  to  note  that  (l)  is  a  very  comprehensive  equation.  It 
covers  the  oomplete  range  of  linear  system  analysis  and  synthesis  covered  by 
Laplace  or  Fourier  transform  methods,  M(t)  may  be  equivalent  to  the  opera¬ 
tion  of  differentiating,  or  integrating,  or,  in  fact,  may  be  equivalent  to 
any  ordinary  integro-differential  equation.  Due  to  this  faot  that  the  super¬ 
position  integral  embraoes  all  of  linear  system  analysis  and  synthesis,  it 
provides  a  proper  basis  for  the  development  of  any  form  of  operational  cal¬ 
culus. 


To  construct  an  operational  oalculus,  one  defines,  corresponding  to  a 
given  time  funotion  F(t),  another  funotion  F(x),  called  the  transform  of 
F(t).  If  then,  corresponding  to  the  superposition  integral  (l),  there 
exists  a  simple  relation  between  the  transforms  of  F(t),  M(t)  and  H(t),  we 
shall  have  an  operational  method.  In  Laplace  transform  theory,  the  trans¬ 
form  relation  is  L(H(t)i  =  L(M(t)J  •  L(F(t)j.  In  Fourier  transform  theory 
the  relation  is  identical  in  form,  being  H(j<o)  =  M(jcc)  F(ja>)«*  In  the  new 
method^to  bo  set  forth  here,  the  transform  relation  is  again  of  this  some 
form,  H(x)  =  ff(x)  •  ^"(x),  where  these  respective  functions  are  the  BPoly- 
nomial-Transforms”  of  H(t),  M(t),  F(t).  We  must  now  define,  of  course,  this 
new  type  of  transform,  which  we  abbreviate  as  "P-Transf  crm1* . 

POLYNOMIAL  TRANSFORMS 


Corresponding  to  a  given  time  funotion  F(t),  we  define  now  a  transform 
T(x),  called  here  a  P- Transform,  by  tbs  equation. 


L. 

F(x)  *P  F(t)  =»2Z»F(nv)xn  (2) 

n=-K 


where  the  F'(nv)  constitutes  a  sequence  of  ordinates  of  F(t)  equally  spaced 
a  time  v  apart.  The  respective  factors  xn  may,  for  the  moment,  be  thought 
of  as  a  means  of  ordering  the  sequence  F(nv)  of  ordinates.  The  word  ordinate 
F(nv)  may,  of  course,  include  an  ordered  sequence  of  empirical  values,  taken 
for  example,  from  a  table.  And  the  independent  variable  t  does  not,  of 
course,  have  to  be  time. 


We  also  define  an  operator  P  "*  ,  inverse  to  P  as  defined  in  (2).  Thus, 
having  given  a  P-tranaform  F(x),  that  is  to  say,  its  sequence  of  ordinates 
F(nv),  we  shall  understand  by  the  expression 


f(t)  =  P"1  F(x) 


(3) 


the  interpolation  of  a  time  function  f(t)  having  the  values  f(nv)=F(nv) 

(n  =*  «K  to  L) »  Graphically,  this  amounts  to  drawing  a  curve  f(t)  exaotly 
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through  the  pointa  F(nv).  We  shall  call  f(t)  in  (^)  a  DETERMINING  FUNCTION 
or  TIME  FUNCTION  corresponding  to  the  P-transform  F(x) .  The  function  f(t) 
ia  an  approximation  to  F(t)  which  can  be  made  exact.  Before  a  determine 
function  f(t)  becomes  uniquely  defined,  we  shall  have  to  prescribe  how  its 
ordinate  sequence  f(nv)  =  F(nv)  is  i-o  be  interpolated.  For  conciseness,  we 
shall  occasionally  writ®  Fn  for  F(nv). 

CARDINAL  INTERPOLATION 

For  pr«»*nt  purpo»*a,  interpolations  will  be  limited  to  the  form 
L 

f(t)  =«  ZZLpnL{t»nf)  (4) 

na-K 

where  the  interpolating  function  L(t)  will  be  taken  as  the  CARDINAL  FUNCTION 

IT  . 


(4) 


L(t)  = 


sin 


(5) 


a  graph  of  whioh  is  shown  in  the  upper  sketoh  of  Figure  3. 
has  the  properties 

,  L(t-nv)  =  1  for  t  =  nvj 

=  0  for  t  =  kv  > 
n,k  =  integers,  n^  k 


This  function 


(6) 


L(t)  =  L(-t) ,  (7) 

is  continuous  for  all  t  and  its  tails  attenuate  hyperbolioally  to  rero  as 
t  approaches  infinity. 

The  meaning  of  the  inverse  transform  f(t)  =  P  F(x)  may  be  seen  dearly 
as  follows.  If  wo  multiply  L(t-nv)  by  I,,  and  graph  this  product  in  proper 
position  on  the  graph  of  F(t),  see  Figure  3»  obtain  a  curve  enclosing  the 
shaded  area.  This  curve  passes  exactly  through  this  partioular  point  (nv,F(1) 
and  oroases  the  time  axis  at  every  other  ordinate  of  the  sequence  .  If  we 
do  the  like  for  every  ordinate  Fn  and  add  the  separate  interpolating  ourves 
so  obtained  for  all  valuer,  of  t ,  we  see  that  the  sum  curve  (I4.)  representing 
the  inverse  transform  f(t)  is  an  approximation  or  interpolation,  continuous 
for  all  t,  of  the  time  function  F(t).  Note  that  this  interpolation  f(t) 
passes  exact ly  through  all  points  of  the  ordinate  sequence  Fu  .  Note  further 
that,  by  passing  to  the  limit  as  v  ^  0,  f(t)  oan  be  made  exactly  equal  to 
F(t)  if  l'(t)  ia  continuous. 
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It  is  also  of  interest  to  note  that  the  area  under  the  curve  FrtL(t-nv) 
is  F0v.  This  is  the  correct  amount  of  impulse  to  be  associated  with  the 
ordinate  Fn  . 

(!>) 

The  cardinal  function  also  has  the  following  useful  properties  .  L(t) 
is  equal,  up  to  a  constant  multiplier  v,  to  its  auto-convolu bion,  namely, 

c  00 

L(s)L(t-s)ds  =  v  L(t)  (8) 

-CO 

From  this  it  is  easy  to  show  that  the  function  S8t  L(t-nv),  where  k,n  are 
integers,  is  orthogonal,  namely, 

r 

L(s-nv)L(  s-kv)  d*  -  v  when  n  =  k,  =0  when  n  j  k,  (9) 

-  OJ 

The  properties  (6)  to  (9)  will  suffioe  for  our  present  purposes. 

PRINCIPAL  THEOREMi  We  are  now  in  position  to  proceed  with  the  statement 
and  proof  of  our  principal  theorem: 


PRINCIPAL  THEOREMi  THE  CONVOLUTION  OR  SUPERPOSITION  INTEGRAL, 


H(t)  = 


f  F(s)M(t-a)ds 


-ao 


IMPLIES  THE  P- TRANSFORM  RELATION 
H(x)  =  v  F(x)  M(x) 


(1) 


(10) 


AND  CONVERSELY,  IN  THE  SENSE  OF  APPROXIMATION  WHICH  CAN  BE  MADE 
EXACT,  THE  TRANSFORM  RELATION  IMPLIES  THE  SUPERPOSITION  INTEGRAL. 


We  now  prove  thia  theorem.  The  ti*»  functions  F(t),  M(t)  and  H(t) 
needed  in  (l)  way  be  represented  by  the  interpolations  of  the  cardinal 
type  (U)» 

F(a)  =»  ZH  F-  L(e-Jv),  F;  «*  0  outside 

5-T  J  J 


U(t-«)  =2^  l(t-i-kv).  * 
H(t)  a  L _ iR  L(t-ny),  H  0 


n~-N 


**  0  outsid*  -X  £,  k ir  L 

outside  -IPb  n  ir  M 


(n) 
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Inserting  thes®  into  (l)  and  commuting  the  order  of  summation  and  integra¬ 
tion. 


H 


H(t) 


=£ 


H  L(t-nv) 


I  L 

Ga  -■* 

Ff-  Z—M 
j»-7  J 1^-K 


r« 


L(s-jv)L(t-8-kv)ds. 


Sotting  x  a  s-Jv  and  invoking  (6),  we  have 

M,  I  l, 

H(t)  a  Z_J>  H  L(t-nv)  =  vZH  F-  X*  M  L  [t-(j+k)v]  . 
n=-H  j  =  ~-j  3  k  : -K  1,4 

Making  uae  of  (6),  the  last  expression  will  yield  explioit  values  for  the 
H0  in  terms  of  the  Fr  and  Thus,  let  us  set  t  =  iv  in  the  last  equa¬ 
tion!  on  the  left,  all  terms  vanish  exoept  In  the  seoond  summation 

on  the  right,  all  terns  vanish  exoept  that  for  whioh  i  =  j  +  k,  the  surviv¬ 
ing  term  being  M.  .  .  Hence, 

t-M  t 

H.  «  v  F.  M-  |  whioh  we  write  as 

Hn”r^FJMn-i  <12> 


This  is  a  very  useful  result. 


To 

for  the 


see  more  clearly  what  (12)  Moans,  let  us  now  form  the  P-transform 
sequent) e  H  i 


aw  . £ H . zi EL*.  • 

n«-H  «  Ow-M  j*-J  i  n'J 


-v’ktjit 

js-T  '  nr-N 


U  .  x 

"'J 


«-a 


Replacing  the  index  n  -  J  by  k  and  changing  the  limits  on  the  n  summation 
to  corresponding  values  of  k,  we  have 

H(x)  a  XI  H  xn  «  v  21  F;  (12'  ) 

n.-M  •»  3trTi  kt^  * 

In  the  applications  of  this  operational  method,  we  shall  frequently  encounter 
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equation#  of  this  form  (12}),  for  which  reason  (12  )  deserves  a  careful 
study. 


In  equations  (11),  we  placed  certain  restrictions  upon  the  F;  ,  Mkand 
Hn,  stating  that  these  ordimto#  vanish  outside  certain  ranges  -  it  amounts 
to  the  same  thing  to  nay  that  outside  these  ranges  those  sets  of  ordinates 
become  negligible,.  For  the  purposes  of  system  analysis  and  synthesis,  w© 
can  place  upon  the  ranges  concerned  the  furtl®r  restriction  that  M  =  I  +  L, 

BJ  “>  <J  +  K.  Ibis  last  restriction  will  not  be  maintained  in  a  later  report 
dealing  with  noise  smoothing  problems,  and  this  will  require  a  more  detailed 
study  of  the  form  (12*) . 

We  can  now  simplify  (121)  *a  follows.  Equating  coefficients  of  like 
powers  of  k  on  both  sides  of  (12* ),  we  have, 

.  M-j  M-I 

“  “  F_  x*  2 _ *  U  2  E  xJ]^L — .  M  * 1  -*•  F  x  Z  K  x*\ 

v  k«?“N<*T  ll  J  Wa-N-i  k  l  b. 


Recalling  the  additional  restriction  that  K  =  I  +  L,  we  s#o  that  the  upper 
limits  on  these  successive  aummationa  decrease  from  M  +  J>LtoM  -  I  z*  Lt 
and  eino®  0  for  all  k  >  L,  all  these  upper  limits  can  be  replaced  by  L» 

In  similar  warmer,  the  successive  lower  limits  go  fro®  -N  +  J  »  -K  to 
-N  ~I<“K,  and  einoo  =»  0  for  all  k<~K,  we  see  that  all  these  lower  limits 
oan  b®  replaced  by  -K«  Hence,  in  (12*)  the  lower  and  upper  limits  k  =  -fl  ~  j 
and  k  =»  M  -  j  com  be  replaced  respectively  by  k  =  -K  and  kn~L,  Accordingly, 
(12* )  become# 


H( 


*)  «  <c__  H  X*  ®  V  F<  * J  iH*  1 
n  j«~T  J  k=-kT 


j«~T 

5=  I  -»•  L  and  N  s  J  t 


Ibis  relation  may  also  be  illustrated  by  multiplying  out  two  polynomials 
of  low  degrees. 


U3) 


We  not®  am  that  the  two  polynomial  f motors  on  the  right  side  of  (13) 
are  respectively  F(x)  and  M(x).  Hence  H(x)  -  v  ?(x)  $!(x).  This  proves  tiie 
first  part  of  our  principal  theorem. 

It  remains  to  prove  the  converse,  namely,  that  (1)  follows  from  (10). 
Lot  us  write  (12)  in  th®  functional  aotationr 

H(iv)  =  2_F(jv)  M(iv-jv)  v 
avst-Xv 

Wow  let  the  time  interval  v  «  ds-wO.  Since  v  is  decreasing,  the  number  of 
intervals  must  increase  to  oover  th*  given  ranges  of  variables*  hence,  we 
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s«t  iva  t,  jv  =  8,  Jv  =  T(  ,  and  Iv  =  T^»  Passing  to  the  limit,  than, 

Iv. 

H(iv)  =  .n  F(jv)M(iv  -  jv)  v 
jVs-Tv 


gives 


/T». 

H(t)  =  |  F(s)  M(  t*»s ) ds o 

T, 


With  proper  restrictions  as  to  convergence  of  this  integral,  already  im¬ 
plied  in  (1),  we  oan  let  T,~v«x>  and  Tz  and  so  obtain  (l)»  H««vn©, 

since  (12)  passes  over  in  the  limit  to  (l),  the  equation  (10)  whim  ^s  equi¬ 
valent  to  (12)  is  an  approximation  to  (l)  which  oan  be  mad®  exact.  This 
oompletes  the  proof  of  our  principal  theorem. 


It  is  interesting  to  compare  (see  Appendix  A)  the  present  operational 
methods  with  the  familiar  methods  based  upon  the  Fourier  and  Laplace  transform 
theory. 


LINEAR  SYSTEM  THEORY 


There  are  three  typical  problems  in  linear  system  theoryi  th®  analysis 
problem,  the  synthesis  or  design  problem,  and  the  instrument  problem. 

In  the  typical  ANALYSIS  PROBLEM,  one  has  the  system  M(t)  given  (the 
memory  function  being  used  here  to  characterize  the  system),  one  has  given 
the  input  F(t),  and  one  has  to  o&loulate  the  system  output  H(t). 

In  the  typical  SYNTHESIS  PROBLEM,  one  has  given  the  input  F(t)  and  the 
output  H(t),  and  one  has  to  calculate  the  system  memory  function  M(t)  - 
aooepting  M(t)  as  the  attorney  or  mathematical  representative  of  th©  system 
itself,  and  leaving  it  up  to  the  design  engineer  actually  to  make  the  system 
so  it  will  have  this  memory  funotio.  response  to  unit  impulse. 

In  the  INSTRUMENT  PROBLEM,  one  has  given  the  system  M(t)  (the  measuring 
device),  one  has  given  its  output  H(t)  (the  measurement),  and  one  has  to  find 
out  what  the  input  F(t)  is  or  was. 

All  these  problems  are  easily  solved  by  the  present  operational  oaloulus, 
nothing  but  the  ordinary  algebraic  operations  of  polynomial  multiplication 
or  division  being  needed. 


Consider  the  analysis  problem.  Having  given  the  input  F(t),  one  picks 
off  its  graph  or  calculates  a  sequenoe  of  ordinates  F»  equally  spaced  at 
integral  time  multiples  jv,  and  fora*  the  input  P- transform 


?(*) 


A5iTR  658I 


7 


In  similar  Banner*  one  forma  (practically  by  inspection)  tha  memory  P-transform 

kx)  »  2-1 . 


On®  than  Multiplies  thee*  two  polynomial a  by  th# 
alg«br®  to  obtain 


H(sc)  s»  ▼  F(x)M(x) 


familiar  algor  it  Jim  of  ordinary 

I-tL. 

r-n 

st  2. - -  H  3tn  « 

<r=-X-K  * 


On®  then  pick*  off  the  coefficients  Hn„  plots  these  at  tl*e  reapeotive  pointe 
nv,  fairs  in  a  tin®  curve  h(t)  tlirough  these  points*  and  accepts  this  ourvo 
as  an  adequate  interpolation  of  the  required  response  H(t), 

In  the  synthesis  problem,  one  forms  as  above  the  P-trans forms  P(x) 
and  H(x)  of  given  input  and  output.  One  then  divides,  by  ordinary  poly¬ 
nomial  division  the  first  by  the  second  to  get  the  memory  transform 


picks  off  and  graphs  the  ordinate  sequence  and  obtains  the  interpolated 
time  function  M(t)» 


In  the  instrument  problem,  one  forms  the  P-transforsas  of  th©  given  mea¬ 
surement  H(t)  and  given  memory  M(t),  divides  to  get  the  polynomial 


<N/ 

F(x) 


1  S2l 

▼  ifc)* 


picks  off  this  P-transform  F(x)  its  ordinate  sequence  Fn ,  graphs  these  or¬ 
dinates,  and  obtains  the  interpolated  measured  quantity  F(t), 


It  will  be  noted  in  all  three  of  these  problems  that,  because  of  the 
obvious  relation  between  P-transform  and  corresponding  time  function,  one 
remains  essentially  in  the  time  domain  itself.  The  great  intuitive  value 
of  this  fact  will  be  made  apparent  in  a  later  paper  (AP  Technical  Report  Ho. 
6586)  in  whioh  the  Kolmogoroff-JViener  theory  of  noise  smoothing  and  predict¬ 
ing  will  be  presented  in  terms  of  the  new  operational  method.  The  essential 
ideas  of  noise  analysis  (regarded  by  engineers  as  too  abstruse  for  practical 
use)  take  on  a  simple,  obvious  character. 


The  value  of  this  new  method  as  a  practical  way  of  getting  numerical 
results  will  be  Illustrated  later  in  a  numerical  example.  This  example  will 
be  delayed  in  order  to  set  forth  a  number  of  useful  theorems. 
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USEFUL  THEOREMS 


There  ara  many  proportie a  of  the  P- transformation  which  are  useful  in 
applications  of  tho  present  operational  method.  In  general ,  these  proposi¬ 
tions  correspond  to  parallel  propositions  in  Laplace  transform  theory.  Some 
aro  sufficiently  obvious  to  bo  stated  without  proof. 


Th*  It  Tho  following  operations,  respectively  in  the  time  and  trans¬ 
form  domains,  correspond! 

If  H(t)  a  F(t)  +  K(t), 

then  H(x)  a  F(x)  +  K(x) 

Th.  2s  If  H(t)  =  o  F(t), 

then,  H(x)  a  o  F(x) 


Th.  ?t  If  fl(t)  a  F(t-anr)  ] 
then  H(x)  -  x'™  F(x) 


(iU) 


(15) 


(16) 


Th©  last  theorem  moans  that  multiplication  of  th©  transform  P(.x)  **  PF(t) 
by  xm  corresponds  in  the  time  domain  to  shifting  the  time  funotion  by 
amount  uw  into  the  future. 


Th.  Us  If  H(t)  =  F(-t) 

then  H(x)  =  2H  f’(nv)  nn 


(17) 


This  means  that,  corresponding  to  the  time  domain  operation  of  reflecting 
the  funotion  F(t)  about  the  time  origin,  on®  merely  has  to  change  the  sign® 
of  the  exponents  of  x  in  tho  transform  F(x)  =  ^7*.  F(nv)xn« 


Th.  5s  If  H(t)  =AF(t)  =  f(t+v)  -F(t)| 
then  H(x)  =  (I  -  l)  F(jc) 


(19) 


which  means  that  the  operation  of  taking  finite  differences  in  tho  time  do¬ 
main  corresponds  to  multiplication  by  Cl  -  1)  in  the  transform  domain. 


Th.  6  s 


If  H(t)  =  t  F(t) 

iv>  dF(x) 

then  H(x)  =  xv  — 


(19) 
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which  moans  that  multiplication  by  t  in  tho  time  domain  corresponds  to  the 
operation  vx  in  the  transform  domain.  If  one  write  a  i(x)  =  F(nv)  x" 

»Ad  differentiates „  the  proposition  readily  followii, 

Th*  7 i  A  theorem  for  integration  in  the  P  domain  goes  aw  followut 


If  H(t) 


*1*1 


t  +  v 


thou  H(x) 


vx 

oo 


r* 

j  F(x) 


dx  y 


whore  F(x)  =  <f — >  F(nv)  xn 
n*o 


(20) 


which  follows  readily  by  substituting  the  third  equation  into  tho  second. 
The  term  F(-v)x~,  would  introduce  a  logarithmic  singularity  and  is  reserved 
for  later  study. 


Tho  amount  of  impulse  H(nv)  duo  to  the  notion  of  a  driving  Amotion 
F(t)  acting  over  a  time  from  0  to  nv  is  the  definite  integral 

vW 

F  ( t )  dt 
'o 

and  is  the  area  under  the  curve  F(t)  from  0  to  t  ~  nv.  Concerning  the 
ooquenoo  of  values  Hn  H(uv)„  we  can  state  the  theorem 


f 


Th.  8»  If  in  the  time  domain 


(  nv 

H(uv)  a  I  F(t)  dt 

•'ft 


in  tho  P- domain 


o 

oo 

V 

itZ  K(x) 

i—X 

V 

2 

“  2 

i*  ..-j 

0  k«i 

03 


where  f 


— 1 

(x)  *2_ 
kao 


K 


(21) 


and  where  H(x)  a  /  H  xn 

'  •  itwww »  . 

IUM3  " 

In  otlier  words  Hr,  is  tho  coefficient  of  x*'  in  the  second  equation  of  (21), 
This  theorem  will  be  useful  later  in  reducing  ordinary  linear  differential 
equations  to  the  polynomial  operations  of  ordinary  algebra.  For  this  rea¬ 
son,  it  deserves  a  careful  discussion, 

Tho  area  under  any  portion  of  the  curve  F(t)  io,  to  a  sufficient  ap~ 

eo^ 

praxlmation,  equal  to  tho  area  under  that  part  of  its  interpolation  / . F  L(t-kv), 
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The  area  from  0  to  t  -  nv  ie  then 

>w 


H(nv)  *  H|(  “  I  r*  v,  L(  t-kv)dt 

*  #  |<ho 


nV 


“  sL.  F  I  L(ti-kv)dt 
V:  ^  o  ** 


For  oonolaonese,  we  write 


o.  «« 
N 


(  L(  t-kv)dt 


uo  that 


..  r..  v  s 


U“  «* 


ki  > 

k  k- 


(iW) 

U'» 


It  remains,  of  course*  to  evaluate  S^» 

If  w<>  \iao  as  interpolating  function  the  triangle  shown  lit  Figure  I;, 
whioh  pr  odiums  the  straight  lino  approximation  or  polygon  shown  in  this  lower 
sketch,  it  is  obvious  (sines  the  area,  of  the  triangle  1«  v)  that 


8^  --i  v  for  Or,  k  *'  a 
■«  ~  for  k  ^  o  or  n 


Note  that  of  the  triangular  interpolating  functions  at  t  »  0  and  t  w  nv,  the 
half  triangles  out  aide  tho  range  0  to  uv  are  excluded  in  order  to  provide 
lie t to v  appr  o.x  i \ua  1 1  on n 

Tim  cardinal  function  (h)  givey  approximately  the  same  result  -  and  thi« 
iti  dismissed  in  Appendix  II  to  avoid  interrupting  tho  continuity  hero u 

Fx pan ding  (11?)  *i'd  usiivg  (I’Ji) , 


K  -  Mo  +  S,  F, 


♦  SJ\  +  f  Vi  Fn-» 


+  B  i' 

M  »1 


2 


F  i  vF  ►  vF 


I  VF 


■*7 


which  Ik  also  obvious  from  a  calculation  of  the  trapocoidal  areas  of  tho 
interpolation  a  a  ehown  in  Figure  1;*  Thie  result  can  be  written  as 


v 

O 


(2F,  *•  8F,  t  2F,  l- 


-  n.)  -*■„] 


(25) 
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Tim  quantity  In  l, ho  pm  rant- he  no  ■  (  )  in  tlui  nod'l'l  o  1  ant  of  a"  In  t)k»  product 

iK> 

of  t.hm  two  polynomial  a  ( J.'.  (1  *•  2  x  i  2x'  i  2x’'  t •••ad.inf «) ,  which 

can  he  readily  mhown  hy  *  bit  of  pencil  work,  Furthermore ,  e»  division  by 

1-x  will  uhow  that.  (1  1  2x  1  2xl  t  :.'x'V  » •••ml. inf.)  -  }-*•  .  The  remaining', 

1>"X 

term  F  in  the  bracket  of  (2b)  in  the  ooef f lolent.  of  x*'  An  the  polynomial 


1.  x  «  1‘ 

h  *»  1 


X  \ 

where  the  term  for  K  «•  0  lu  omitted  to  ft  void  cancel  lint'. 

*  1-x 

the  x°  term  of  t Ho  previous  perl.  Hone*,  it  appears  that  11  „  la  the  coef¬ 
ficient  of  x"  ii\  the  1*~  transform. 

IKl  V*’ 

.  v  lu  V  k  K  V  \  *  Iv 
HU)  -  -*-■  /  .  xK 


V  1  <  X 
iff  T-x 


I-  1 

wo 

F{x)  -  *  K„  >.*  xk 
•  k*i 


v  1  i  x  7; ,  \ 

“?  \~i  1,(10 


K,  vx 


(in*) 


from  which  Theorem  M  full own* 

A  better  approximation  than  (2b)  to  the  xreu  or  impulse  H<(  1b  given 
by  Simpson' ■  one-third  rule 


(hi 


“n  “ 


3 


[Va  *  ]»K.  *  VV>.  (  J*Fw  *  iM%  +  ,Vm  .1 


0*6) 


where  n  i<\  even,.  By  direct  execution  of  the  opera  Ilona  below,  it  can  be 
nhown  that 

x.  o'* 

Of,  v  V  iWjxtx  7,,  v  V  V"  h 

su) n»)  -  -  >^;r  *  • 

k  mvii.u 

Hence,  ( 2? )  may  he  need  in  place  of  the  second  equation  of  (21)  -  one  ahould 
recall  then  tluat  a  is  to  be  taken  only  uh  an  even  number* 


0*7) 


A  large  olano  of  interpolating  funotioiw  from  which  additional  approxi- 
wationw  may  lie  derived  aro  given  by  Schoenberg  64), 

If  one  write e  the  first  equation  of  (21)  in  the  form  H(t)  a  |^K(t)dt, 

it  followa  that  K ( fc )  -  H(t)  = In  thiu  oaae,  the  oonatant  lower  limit 

plays  no  part,  According,  we  may  drop  the  term  in  t’0  fr  om  the  second  equation 

nrt 

of  (21),  Solving  thin  oeoond  equation  for  H(x)  =  F(x),  wo  have  the  theorem. 
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Th.  9 1  Corr«8j)omllnK  to  the  time  domain  operation  of  ditToroni.lui.ion, 


ft<t)  -  ^ 


dt 

<r,  j>  ,•<  i 

we  Iwvo  in  the  I’-domaln  H(x)  ^  n(x)  f 

where  !!(«)«?  .  H(kv)x^ 

K“*> 


A  nt.  ill  better  l’-domain  operation  corresponding  to  differential.!  on  in  tho 
tilin'  domain  in  obtained  by  a  Mitnilar  inversion  ol‘  ( ) ,  thus 


5  1-x 

4t  „ 

v  it!  i 


JX»X  1 


ti(«) 


and  thin  may  bo  used  in  pin chi  of  the  second  equation  in  («.Hl)« 

Tho  no  theorems  might  bo  continued,  mono  or  1  o  ;i  «  paral.lol.lnp,  tho** 
of  Laplace  transform  theory.  Knough  have  boon  given,  however,  to  sketch 
out  the  procedure  for  developing  tho  now  operational  calculus. 

An  example  (coming  under  Th.  II)  will  now  bo  worked  out  to  illustrate 
the  arithmetic  convenience  of  the  new  methods* 


NUMLK1CAI.  fc.XAfcU’Iii 


Ab  an  example  under  Th.  U,  wo  calculate  the  ordinates  H(nv)  of  th* 
nv  v 

integral  sin  t  dt.  Taking;  the  time  interval  us  v  ™  a  18°,  the  P~ 

■'ll 

transform  of  sin  t,  taken  from  a  trigonometric  table,  is 

F(x)  ”>  0  +  «,31x  +  i  ,Blx'  t-  «95xf  +  x  f  ,95x  +  etc* 

From  (21),  tho  transform  operator  i«  53  *157  y“'x  0  Hence,  noting 

that  1  (0)  »  ulri  0  a  0,  wo  have 

H(x)  a  0,137  ““  (0  t  .jlx  +  859x"'  +  .eix4  +  .95*'  +  lx4  t  .95*  U  +•***) 


With  a  bit  of  pencil  work  on  tho  side  (one  may  multiply  fir at  by  1  +  x  and 
then  divido  by  1  -  x  or  vise  vorea),  tho  reader  will  readily  verify  the 
results  calculated  hero,  the  second  gives  valuew  calculated  from  the  analy¬ 
tical  solution  ,,w 

I  «in  t  dt  •-  l-oos  nr. 
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H  v 

H, 

H* 

Hr,-  H, 

Hy 

H,a 

Ha 

0 

<•05 

.19 

Jll 

.69 

.98  1.30  1.57 

1.79 

l*9h 

1.98 

1.9U 

1.79 

0 

.03 

.19 

.m 

.69 

1.00  1.31  1.59 

1.61 

1.93 

2.00 

1.95 

1.81 

These  results  aro  nuff ioiontly  accurate  to  illuatrate  the  simplicity  and  pre¬ 
cision  of  the  new  methods.  A#  an  additional  bit  of  prnotioe,  the  reader 
may  use,  instead  of  the  integrating  operator  above,  the  bettor  one  given  in 
(27),  namely. 


v  l+Lpc-t-x*' 
3  1-x  ' 


.101+7 


l-flpc-fxv 

l-xK 


a 


A  TABLE  OF  P- TRANSFORMS 

For  renaonn  brought  out  00  clearly  by  Gardner  and  Barnes*  in  the  oa«« 
of  Laplace  trami form  theory,  we  find  here  that  a  table  of  P-tranaf orjnw  ar® 
useful.  A  few  oxamplen  only  can  be  worked  out  her®,  and  these  will  bo 
limited  to  the  elementary  function*  usually  encountered  in  linear  eyatete 
theory.  The  oaloulationa  amount,  of  couroe,  to  diroot  application#  of  the 
definition 

td 

P(x)  a  P  [P(t)]  -  F(nv)  xn  (2) 

VS=:~N 


of  the  P-trnnafom. 

The  unit  atop  function  S(t)  is  defined  by  the  equation 

S(t)  a  1  for  1 3?  0 
0  for  0 

The  corresponding  P-transform  S(x)  iss  S(.x)  a  1  +•  x  +  xx  +*«.+xn  +  etc 

"=  2—<.  xn  ,  But  1  +  x  +  x,1'  +•  ..adeinf  o  =*  ,  henoo 

n«o  1*»JK 


in  tk  closed  and  simple  form. 

The  truncated  exponential  function  is  defined  by 


(30) 


(31) 


AFTR  658I 


E(t)  =  e  ^  for  0 


=  0 


(32) 


<0 


—  anv  „o(v  \H  ,  k  — >  —otW  r> 

The  ordinate  sequence  ia  E(nv)  =  e  =  (©  )  j  hene®,  E(x)  «  /  ®  x 

OO 

sr  ,  .  n 

=  / _ „  (©  jc)  *  Accordingly,  we  oaa  writ® 


n~a 


not* 


oo 


S(x)  =  Z1  (axf  =“5~" 

n-o  *  ** 


(33) 


where  a  =  a 


v 


Wot©  that  by  setting  aL  =  0,  and  hence  a  »  1#  we  recover  th®  previous  re- 
suit  for  th®  unit  step  function, 

aft  *co,t 

In  similar  fashion,  by  calculating  P  [  e™  ©^  land  asperating  real 
and  imaginary  parts,  we  find  for  the  truncated  damped  sinusoids  the  poly¬ 
nomial  transforms „ 


r  —  ot  fc 

P^o  sin  tut 


1  U 

i  T-Y2a 


(a  sinoiv)  x 


(2a  oo*unv)  x-Hf'x1' 


m 


where  a  =  e 
...  alt 


•  civ 


„r  .  ~l  l-(a  ooncoy)  x 

P[  e  oos  u,t  I  *X^STooSS¥T^x^ 


(33) 


Application  of  Theorem  6  (Equations  19)  will  introduce  a  factor  t  into 
any  one  of  the  above  functions,  corresponding  to  which  one  applies  to  the 

P~t ran* form  th©  operation  vx  4~  •  As  an  example  of  this,  consider  the  ramp 
function 

R(t)  «  t  ?(t)  (36) 

Its  transform  ia  R(x)  =  vx^~-  S(x)  -  v/-^L.  (t^t)*  h®bo® 


ar  h-5 


R(x)  a  Y’jT'^yx,3  ^(C+x+2xXt«*«+nxn  +  etc) 

Tli®  respective  coefficients  0,  v,  2v#0..,nv,  etc,  in  the  last  expression 
are  obviously  the  sequenoe  of  ordinates  of  this  ramp  function. 


(37) 


IsSultipl©  application  of  Theorems  6,  1  and  2  to  the  above  functions ‘and 
their  transforms  will  yield  the  transforms  of  any  polynomial  in  t  times  the 
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3  x.  -  5T 

original  funotioa  -  suoh  as  (2t  +5^  -3)  ®  •  This  corresponds,  of  course, 

to  a  multiple  pole  (degenerate  eigen  functions)  on  the  complex  plane  of  the 
Lapl&oe  transform. 

The  unit  impulse  S(t)  at  time  t  =  0  is  defined  by  the  conditions 
S  (t)  =0  for  t^O 

r*>  (38) 

6  (t)  dt  --  1 

*Lo 

CO 

Prom  (1+)  we  can  write  (t)  =  2  ^  <S  (kv)L(t-kv)  which  (since  6  (kv)  =  0 

!<•-•■  -<*> 

for  0)  reduces  to  S(t)  =  $  (0)  L(t).  Using  this  result  in  the  neoond 

,  CD 

S  (t)  dt  =  8  (o)  I  L(t)dt  =  vS(o)  =  1*  Hence 

1  -»>  co,  , 

6(0)  =  — ,  Insertion  of  this  into  the  transform  S(x)  =  22  S  (kv)  xK  give* 

▼  -co 

S(x)  =  1  (39) 

The  time  origin  of  any  of  the  above  functions  may  be  changed  from  0  to 
mv  by  use  of  our  "shifting"  Theorem  3»  The  corresponding  operation  in  the 
P-domain  ic  multiplication  by  xm ,  For  example,  the  transform  of  the  delta 
function  located  "at”  mv,  namely,  S  (t-mv)  is 

n  xm 

P  [  6 (t-mv) j  =  ~  (39  ) 

These  results  may  be  oolleoted  in  a  table,  one  column  listing  the  time 
function,  the  ether  the  corresponding  P-transform.  Suoh  a  table  can  b®  made 
quite  extensive  and  would  serve  the  same  general  purposes  as  a  table  of 
Laplace  transforms  or  a  table  of  Fourier  transforms.  The  few  items  listed 
above  suffice  for  our  present  purpose  of  showing  how  suoh  a  table  can  be 
constructed. 

Wo  now  apply  these  methods  to  the  solution  of  ordinary  linear  differen¬ 
tial  equations,  , 


equation  of  (38)  gives  J 


LINEAR  DIFFERENTIAL  EQUATIONS 

The  few  theorems  set  down  above  enable  us  to  solve  ordinary  linear 
differential  equations (  The  general  method  of  att&ok  can  be  olearly  indioated 
by  a  couple  of  simple  examples. 

Consider  the  first  order  equation  with  its  initial  condition. 
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■gf  +  BH  =  F(t) 
with  H(0)  =  H0 


(Uo) 


Integrating  from  0  to  t, 


L 


H(t)dt  +  B  H(t)dt 


(  F(t)dt 

*  O 

=  H(t)  -  H(0)  +  B  J  H(t)dt 


Inserting  the  initial  condition  H(0)  =  Hc  and  write  this  as  H0S(t)  where 
S(t)  ia  the  unit  step  function,  we  have 


t 

H(t )  +  B  f  H(t)dt  = 


J  F(t)dt  +  KoS(t) 

O 


m 


which  is  now  in  proper  condition  to  apply  the  P-transf ormation. 

Theorems  1  and  2  show  that  PC  ]  i*  a  linear  operation  (applies  distri- 
butively  and  commutes  with  a  oonstant) j  henoe,  we  have 

P  [H(t)]  +  BP  [  f  H(t)dt]  »  P  [  {  F(t)dt]  +  H,P  [S(t)]  (1*2) 

**  J  9 

The  first  term  on  the  left  is  simply  H(x).  The  second  term  on  the  left,  by 
Theorem  8^  is  b[~  H(x)  -  *  ^e  last  term  of  which  is  a  minor 

correction  (vanishing  in  the  limit  as  v~**0)  and  will  be  omitted.  The  first 
term  on  the  right  is  similarly  ^  F(x),  The  last  term  on  the  right, 

rw  H- 

see  (Jl),  ia  H  S(x)  =- —  »  Putting  these  values  in  place  of  th®  respective 

/NJ 

tome  in  (Li2),  and  solving  for  H(x) ,  we  have 


3(x) 


2 


1  + 


1+x 

Toe 

FTT 

T  Toe 


Fix) 


+ 


Ho 

jZx 


JTx 
1— x 


(U3) 
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The  coefficient  of  F(x)  corresponds  to  the  transfer  function  of  Laplaoe 
transform  theory  and  to  th®  frequency  response  of  Fourier  transform  theory. 
The  second  tern  in  (1*3)  is  due  to  the  initial  condition.  In  the  so-called 
normal  solution  (initially  quiescent  system  with  dependent  variable  and  all 
its  derivatives  *ero)  the  second  term  in  (1*3)  vanishes  -  whioh  is  usually  the 
case  treated  in  Laplaoe  transfos*m  theory. 

It  is  interesting  to  note  the  correspondence  here  with  usual  operational 
methods.  Thus,  writing  (1*0)  in  operational  notation,  pH  +  BH  =  F,  and  solving 
in  the  form 


H  = 


we  see  that  our  "transfer”  operator  in  (1*3)  could  have  been  obtained  simply 
by  regarding  A  as  our  integrating  operator  —  or  p=~  as 

our  differentiating  operator. 


But  equation  (1*3)  oan  be  simplified  still  more.  Clearing  out  tha  1-x  and 
simplifying,  we  have 


H(x) 


;(.i+»)f(»)+  a?- 

(Bvt-2)+(bv-2)r. 


(U3  ) 


which  is  an  easily  manipulated  form. 

At;  this  stage  in  our  process  of  finding  H(t),  one  could  throw  H(x)  into 
a  sum  of  elementary  types  to  be  found  in  a  table  of  P-traneforma,  and  writ® 
out  the  inverse  transformations  to  obtain  H(t).  The  procedure  here  would 
follow  the  familiar  pattern  of  Laplaoe  transform  theory.  But  this  procedure 
is  not  necessary. 


An  outstanding  virtue  of  this  P-transform  method  lies  in  its  easy  handling 
of  numerical  problems.  For  example,  in  the  oase  of  th®  above  differential 
equation,  every  quantity  except  x  on  the  right  side  of  (1*3  )  would  be  given. 

The  denominator  in  this  expression  >  an  be  (by  dividing  by  Bv+2)  put  into  the 
form  1+Ax.  The  numerator  reduces,  after  collecting  coefficients  of  like  powers 
of  x,  to  a  simple  polynomial.  After  division  by  1+Ax,  we  get  then  a  simple 
polynomial  H(x)  =*  H,  +  H(x  +  H,x*“  +***  +  Hnxn  where  *11  the  Hn  are 

saere  numbers.  This  seqpeno©  .  H„  H , ,  H»,  .  ..,  Hn  ,  , »  •eto.,  of  coefficients  giv® 
tha  suooeasive  ordinates  of  H(nv)  =»  H(t),  and  may  be  plotted  to  giv®  a  graph 
olosely  approximating  our  desired  H(t).  Thus,  no  formal  work  whatever  is 
required  to  pass  from  H(x)  to  H(t)  -  a  more  inspection  of  ft(x)  gives  H(t), 
and  vice  versa. 
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T}*>  above  procedure  for  solving  a  linear  differential  equation  subject 
to  its  initial  conditions  may  be  formalised  into  the  following  steps* 

(1)  Apply  to  the  differential  equation  the  definite  integration 

(  )dt,  as  many  times  as  the  order  of  the  equation.  This  intrcduoes  the 

Jo 

initial  conditions  into  our  solution  and  puts  the  equation  into  proper  form 
for  applying  ova*  P- transformation. 

(2)  Apply  to  the  result  of  Step  1  the  P-transf ormation,  sake  use  of 
Theorems  1  to  8  to  reduce  every  term  containing  H(t),  its  derivatives  or  inte¬ 


grals  into  terms  containing  H(x). 

(3)  Apply  to  the  result  of  Step  2  the  indicated  operations  of  ordinary 
algebra  including  polynomial  multiplication  and  division,  thus  reducing  H(x) 
to  the  simple  polynomial  form  H(x)  =  H0  +  H(  x  +  Hjx3*  +  etc. 

(U)  Pi.k  off  the  coefficients  (or  time  ordinates)  Ho,  H( ,  HI#  H3,  etc., 
plot  at  the  successive  points  0,  v,  2v,  3v,  eto.  on  a  time  axis,  fair  &  ourve 
through  these  points  to  obtain  an  approximation  to  the  required  solution  H(t). 

Consider  now  the  second  order  equation  with  normal  boundary  conditions 
a 


BH  =  F(t) 


m 


with  ll0  ~  H0  =  0 


The  system  is  assumed  to  be  initially  quiaaoent  (dependent  variable  and  all 
its  derivatives  being  zero  at  t=0)  to  conform  to  the  usual  oaae  treated  in 
Laplaoe  transform  theory.  All  this  assumption  does  is  to  omit  a  number  of 
additional  terms  T-hioh  can  be  easily  handled  but  which  clutter  up  the  equations. 


Applying  step  1,  we 


integrate  the  equation  from  0  to  t,  with  the  result 


f  F  ( t )  dt 
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where  w*  omit  «  term  (H„  +  AH0)8(t)  from  the  right  side  due  to  initial  values, 
because  *»  H  «  0*  Integrating  a  second  time. 


x  r*  ra  rc  r 

H(t)  +  A  j"  H(t)dt  +  Ey3rjM(*)d»  =  j  dr  J  F(a 


(Ut4) 


idler*  ft  term  H  fl(t)  la  again  omitted  because  H0  vanishes. 


As  in  step  2  above,  wo  now  apply  the  P- transformation.  The  only  new 
item  her*  l*  the  double  integral.  This  ia  handled  by  a  double  application  of 
Theorem  0.  Thus,  writing  out  the  result  in  full,  we  hav® 

.,t  y>  /U  -r  r\  O 

r  I  I  .  .  v  1+x  f  v  1+x  N  v  F„x  vx  [ 

p[  dr  *g  r_~-  [  g  i~Hx)  -2  i~x  -^0=57  L  (  } 

o  ^  JO 

Hi®  last  tom  van! ah* si  (unless  wo  have  an  impulse  at.  t«0,  a  case  we  omit 
h*r*  for  brevity).  The  last  term  inside  the  braoket  is  only  a  small  correc¬ 
tion  (vanishing  as  v*0)  which  may  be  carried  along  by  thooo  who  wish,  but 
whloh  will  be  omitted  here  for  brevity.  Thus,  wo  shall  have 


p[  j  dr  j  F  ( a )  d  «  j  =  (ey  F(x) 


A  still  better  transform  for  a  double  integral  is  obtained  by  using  the 
transform  operator  of  (2?) .  Let  ua  compromise  here  and  us®  both,,  th®  operator 
of  (27)  once  and  that  of  (21)  once.  This  gives 

p[  J  dr |  F(e)daj  «  y 


V  "  HilpC'bX 


wii«retf  as  above,  we  hav*  owitted  one  term  which  vanishes  and  another  which,, 
for  sufficiently  small  v,  is  negligible. 

Putting  the  result  (U5)  in  for  double  integrals  „  as  well  as  the  more 
familiar  transforms  for  tha  otiwr  terms,  equation  (1x4  )  transforms  into, 

S'(.)  +  A  (£  .£*-)  H(x)  +  B  (J  f  H(x)  =  (’  f  r(x)  . 


Solving  this  for  Il(x)  gives 
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H(x)  = 

1  +  + 


The  parallel  hero  between  the  polynomial  coefficient  of  F(x)  and  the  corres¬ 
ponding  transfer  function  of  Laplace  theory  will  be  immediately  recognised* 

v  l+x 

In  place  of  the  operator  2  Ui3®^  (^6)#  one  IH®y  use  the  operator 


v  1+lix+x 


of  equation  (27)  -  or  om  my  use  ^  &o  tits  first  order  opera¬ 
tor  and  the  operator  of  (U5  )  for  the  second  order# 

By  cle&riag  the  denominators  2(l~x)  from  (I46)  „  we  immediately  simplify 
this  equation  to 


v  (l+2x-<-x  )  F(x)  _ 

H^x'  =  x  +  (I^-2Av+Bvi)  x'x 


(U6') 


Further  reduction  of  H(x)  to  the  simple  polynomial  form  H(x)  =  Hn  +  H,jt  +  HAx'+  eto. 
by  polynomial  division  would  appear  to  be  profitable  if  we  have  numerical  values 
to  put  into  (146').  The  time  function  H(t)  is  immediately  specified,  as  before, 
by  its  ordinate  sequence  Hu,  H(  ,  H?  ,  etc*  picked  off  by  inspection  of  1?(x)« 

The  P-transf orm  solution  of  the  two  simple  equations  worked  out  above 
was  given  in  considerable  detail*  These  examples  should  make  the  pattern 
clear  so  that  one  should  be  able  to  write  down  by  inspection  the  transform 
equation  for  the  normal  solution  of  a  linear  equation  of  any  order*  As  a 
final  example,  we  write  out  the  result  by  inspection  for  the  third  order 
equation 

d^H  ,  d  H  R  dH  PH  _  *,/,  \ 

Tt*  A  It7  +  B  dt  *  CH  ~  f 


or  in  operational  form 


(p J  +  Ap  '  +  Bp  -4-  0)  H=F  • 


From  this. 


p  '  +  Ap *  ■  +  Bp  +  0 


in  which  we  can  now  replace  p  by  our  differentiator  ~  or  one  of  its 

v  1+x 

substitutes*  One  may  also  threw  the  equation  into  the  form 
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H  = 


<V  F 

V 


1  +  a(2)  +  B(I)X  +  o(i)3 
P  P  P 


1  vr  |JuV 

and  replace  by  the  integrator  ^ 


or  one  of  its  substitute 


8* 


In  case  one  is  interested  in  a  problem  where  the  initial  values  of  H 
and  all  its  derivatives  do  not  vanish,  the  more  careful  procedure  used  in 
solving  our  first  equation  above  should  be  followed  in  order  to  get  those 
initial  conditions  introduced  properly  into  the  problem.  Additional  oar© 
must  be  taken  if  one  wishes  to  carry  along  the  small  oorreotions  which 
wore  dropped. 

In  order  to  emphasise  again  the  practical  arithmetic  advantages  of  our 
P-transform  method,  a  numerioal  example  is  worked  out  below. 

EXAMPLES  The  following  example  illustrates  the  arithmetic  oonvenianoo 
of  the  polynomial- transf orsn  method  -  us  well  as  the  other  advantages  in¬ 
herent  in  all  operational  methods. 

ce) 

Consider  the  ^simple  sorvomaehanism  shown  in  Figure  5.  The  torque 
equation  io  T  =  J  sj,  +  F0O,  the  controller  equation  is  T  »  KB,  and  the  error 
equation  is  E  -  9{  -  0^*  Eliminating  0,  and  T,  w®  obtain  for  the  error  the 
differential  equation 

JE  +  FE  +  KEsj2£  +  F&{ 

Transforming  parameters  in  the  usual,  way  by  writing  ujm  =*Vk/J  «  circular 
frequency  of  free  undamped  vibration,  Fc  =  21/  KJ  **  critical  damping  coef¬ 
ficient,  and  C  =  F/Fc  =  damping  ratio  (which  for  conciseness  we  take  here 
as  1),  the  differential  equation  becomes 

E  +  2u^E  +  E  a*  0^  +  2  0' 

The  operational  form  of  this  equation  io 


E  a  -JLZJ&SL--.  e. 

+  2WV  P  t"  ^ 

In  order  to  moke  the  problem  numerioal,  let  us  take  =»  20,  and  the 
driving  funotion  as  the  step  volooity  0^  »  u~>i  tS(t)  =  tS(t).  Then 


p1-*  i^Op  t  ij.00 


(U7) 
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Th*  analytical  solution  of  this  aquation  (normal  solution  with  K(0)  hi  k(0)  »»  0) 

i* 


E 


-  a  (i  +  lot)  ] 


(10) 


2t  is  thin  aquation  (ij?)  which  we  solve  here  by  the  P- trail*  form  method.  1 V« 
■hall  then  compare  our  numerical  result*  with  thorns  of  the  analytical  nolution 
(10) . 


take  the  time  interval  an  v  «  0.05.  Our  differentiating  P- ope rat or 
F- operator  boo  omen  ~  “  )*0  .  Th*  ordinate  sequence  of  the  driving 

l'unotion  le  t-»*0.05  (0,l,2,3,..#eto.) «  henoe,  the  transform  of  the  driving 
i  unction  in  (x)  »  0.05  (^CH-x+J?x  +}x  +I4X  +.«„oto»  ]  ^  P[t].  Inserting  these 
into  (I4 7)  givw* 


0 « 05  ( *+2x  ■*•3**  +  • » •  e 1 0 1> ) 


E3 


1  (l-x)X  +  (l*x)  (1-x) 

5 


(x+2*1-  ■f’lx**. .  .etc , ) 


jjj  q  l  ( l“*x )  (x~f2x  t » «otc » ) 

4  9  -  6ot '+  *'■*■ 

3  v 

.  ,  x  +  x  +  x  +*.eto, 

3  °0  — T^rnrr-- 


-  .010*  t  .07^+  .069X1  +  ,096x^+  ,096xr+  .099x^  UOOx4® 


Tli®  coefficients  hero  arc  the  respective  ordinates  of  the  required  time  func¬ 
tion  E(t)  at  »05,  .10,  .15,  .20, ... eta,  Hi©  following  table,  which  compares 
those  reaulta  with  those  calculated  from  the  analytical  solution  (Ij.8),  shows 
the  practical  degree  of  aoouraoy  of  our  polynomial  transform  method. 


TIME 

O 

O 

• 

.05 

.10 

.15 

.20 

.25 

•  JO  • • • • » 

<ao 

ANALYSIS 

,000 

.oia 

0 

0 

'04 

.087 

.094 

.097 

.099  * • . 0 

.100 

P- TRANSFORM 

•000 

.0)44 

.07U 

.089 

.096 

.098 

•  009  • *  * • 

.100 
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Figure  (>  »lu>w«  h  graph  or  the  driving  function  <bj(t)  i,  un\l  shown 
how  tinny  Vt.  Iti  to  write  ilnviu  by  Inspection  Mu*  driving  I. I’ntwi form  t.'|(x)  *- 
«0<ix  i  .  1  Ox  **  (  ,  Ihx^  >  et.o  .  <-=*  .0*.  (x  i  i  lx* i.. ate,  Thin  figure  nltto 
shows  a  graph  of  the  analyt lonl  solution  together  with  point,*  showing  the 
ordinate*  calculated  above  by  tlw*  F-traiuform  method. 

Tlw*  above*  dissertation  on  linear  differential  equations  ought  ( o  he 
giatttly  el uhorated «  Our  short  treatment  nerves,  however,  t,o  null  Ine  how 
t.ho  new  calculus  may  bo  applied  l,o  linear  dl  fferent  Iti  1  «*pitt  t.i  ons ,  and  out¬ 
lines  how  a  mot* o  complete  preaant.at Ion  o«n  be  worked  out.* 


CONCLUSION 

The  polynomial  -  transform  method  outlined  in  this  paper  iu  worthy  oi* 
a  complete  development,  It  iius  U»  virtue  of  great  simplicity.  The  onulyuls 
follows  thei  imtural  moil®  of  thought (  since  it  remains  always  in  the  time 
domain  (the  l’-traim  f  orm  itself  bo i ng  so  obviously  connected  with  the  time 
domain).  The  very  real  difficulties  in  Laplace  And  Fourier  transform  theory 
in  oalaula ting  dirent  and  Inverse  transforms  tuxi  Avoided  completely  here,  the 
connection  between  time  function  and  lie  P- transform  being  so  obvious  that 
If  one  hu«  either  function  one  can  write  down  the  other  by  inaj.iooti.on.  The 
arithmetic  oonvenionoo  of  the  new  methods  is  obvious,  since  actual  calcula¬ 
tions  involve  only  a  high  mo ho ol  knowledge  of  algebra.  And  finally,  whereas 
In  the  older  methods  the  synthesis  end  instrument  typos  of  problem  are  corn- 
monly  regarded  as  much  more  difficult  than  the  analysis!  problem,  a  note¬ 
worthy  feature  of  the  polynomial- transform  method  in  the  trivial  distinction 
between  these  typos. 

A  paper  will  appear  lator  applying  the  new  methods  to  non-linear  problems. 
Useful  work  on  non-linoar  problems  has  already  been  done  by  Madwod^. 

A  paper  will  appear  shortly  applying  the  P~t  rami  fora  nusthodn  to  "noiso" 
problems  (random  disturbances  in  electrical  circuits,  sorvomeohanioms ,  ill— 
Btsnuuflnts,  and  so  on).  The  Wiener  -  Kolrnogoroff  theory  of  smoothing  and  pre¬ 
dicting  filters  will  bo  reduced  to  a  simple,  immediately  obvious  set  of 
Algebraic  operations. 
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Consider  U\»<  l‘~  transform 


hi 

F  ( x )  “J  /  *  F(nv)  x" 
n«  N 


Lot  nv  »  t,  v  itt  ♦  (),  Nv  *■*>  ,  and  sot 

i.  Oiw 

V  -•*  «  , 


then,  passing  to  th#  limit. 


^  VoUv'  rt  icoA'ov) 

vF ( x )  a  v/-~.vF(nv)  (  «?  )  “  /  <  F(nv)  fi 

"H  -Mv 

(°°  iuit  A 

\  F(t)«S  dt  -  F(oJ) 


•  v 


“OO 


whore  F(«i)  ie  the  Fourier  transform  of  F( t) .  Thus,  we  so®  a  0I030  rolatlon 
between  our  P~transf orm  and  tSie  Fourier  transform.  Had  we  taken  x  =»  e~sv 
and  kept  the  summation  unilateral  from  0  to  to  ,  the  limiting  fora  would 
have  boon  the  Laplace  transform. 


Consider  now  the  inverse  Fourier  transform. 


F(t) 


JL„ 

2Tr 


ao 


A  -lu>i 

F(u))e  dto 


-OO 


This  is,  &s  shown  above,  the  limiting  value  of 


f(t)  =  JL 

2  IT 

ice  v 

as  v-^0*  Replacing  x  bye 
integration,  we  have 


,  4r -  /  v  n  x  •lost 

(v^L.  F(nv)  x  d t-o 

-H 


V 

and  interchanging  order  of  summation  and 
E* 

H  C  y  ) 

XI  F(av)  I  e  d<xi 

— M  »LtT 


V 


V 


V 

2ir 

N 


£  ..  k  F(  rtv ) 


<? 


i(mr-t) 


V“"> 


A-  -s  F(nv) 

n>s-N  —  (t-nv) 

v 


whioh  is  a  determining  or  liiuo  Amotion  corresponding  to  the  P»trenwfom 
^(x).  Thu  a,  tha  invar ae  Fourier  transform  leads  to  tlw*  inverse  P-tranafom® 
Conversely,  one  om  atari  with  tha  inverse  P-tr&nafora  and,  upon  pasoing  to 
tha  limit,  obtain  tha  inverse  Fourier  transform* 
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APPENDIX  li 


In  tiie  proof  of  Theorem  8 ^  the  cardinal  function  itself  raay  bo  used 
instead  of  the  triangular  interpolation  function.  Thus,  in  calculating  the 
coefficients  S^,  of  equation  (22),  we  act 


S 


k 


,.nv  -r 
I  a  in  -•  (t-kv) 

-j;  (t-kv) 


at 


This  can  be  written  as 


„  _  v 


-t.»o 


ain^™  (t-kv) 
Ij;  (t-kv) 


D  [  ~  (t-kv)] 


Changing  the  variable  of  integration  to  u  (t-kv)  and  replacing  the 
limits  by  the  corresponding  values  of  u,  we  have 


S  =2 

k  rr 


-lor  “  ^  L  J- 


ain  u 
u 


du  -J- 


Air 

J  ^  *] 

** 


=  ff'  {  S1  t  ( n-acj'rrj  +  si  (kiv)  | 

rK 

The  a ine-i integral  function  (Si  x  =  ~~ -  du)  increases^  from  0  at 

Jo 

x  =  0  and  oscillates  with  decreasing  amplitude  about  the  value  vy  ,  approach¬ 
'll' 

ing  5  a»i  increases,  the  greatest  overshoot  being  about  c#  with  the  am¬ 
plitude  of  overahoot  decreasing  approximately  hyperbolloally. 


With  this  information  about  the  sine— integral  function,  consider  then 
the  result, 

~  j  Si[(n-k)ir]  +  si(kir)  j-. 

For  k  =  0,  the  seoond  term  vanishes  and  S0  ~~  Si(uTr)  which  approxis*t.,8 
v  11'  v 

if  *  2  =  t  ®vem  for  moderate  values  of  n.  Similarly,  for  k  =  n,  the  first 
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yf 

term  van l oh a a  and  “  9v‘  whioh  oan  again  he  taken  as  y  .  For 

k  «  •$  ,  n>  Sl(^ir)  ,  «c  that  S,^  >-»  v  approximately.  In  fact,  duo  to 
tho  rapid  lnornaao  of  Si(x)  from  0  to  ^  »a  x  inoroaaea  from  0,  all  In¬ 
termediate  values  of  oan  bo  taken  aa  approximately  v.  Accordingly,  we 
liave  appr ox i  umbel  y „ 

S.  =*  v  for  0<  k  <  n 

K 

-  ,Z  for  k  =  0  or  n 

confirming  our  result  (2i*) . 
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FIGURE  I 

Black-box  Representation  of  a 
Physical  System 
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FIGURE  6 


Servomechanism  Input  and  Error  Spectra 
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